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The linear cosmological perturbation theory of an almost homogeneous and isotropic perfect fluid 
universe is reconsidered and formally simplified by introducing new covariant and gauge-invariant 
variables with physical interpretations on hypersurfaces of constant expansion, constant curvature 
or constant energy density. The existence of conserved perturbation quantities on scales larger than 
the Hubble scale is discussed. The quantity which is conserved on large scales in a flat background 
universe may be expressed in terms of the fractional, spatial gradient of the energy density on con- 
stant expansion hypersurfaces or, alternatively, with the help of expansion or curvature perturbation 
variables on hypersurfaces of constant energy density. For nonvanishing background curvature the 
perturbation dynamics is most suitably described in terms of energy density perturbations on hy- 
persurfaces of constant curvature. 
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CS| ' I. INTRODUCTION 



Understanding the dynamics of small inhomogeneities in the early universe is essential for any theory of cosmological 
■ structure formation. The conventional, linear cosmological perturbation theory relies on a splitting of all metric and 
| matter quantities into a homogeneous and isotropic zeroth order and small, first-order perturbations about this 
background. This procedure implies a splitting of the spacetime itself. The nonuniqueness of this splitting is the 
source of the so-called gauge problem A perturbation analysis may either be performed in a specific gauge, the 

traditionally favoured gauge used to be the synchronous gauge (see, e.g., ||), or in a gauge- invariant manner. While 
the choice of a specific gauge corresponds to a specific identification of the points of the fictitious background spacetime 
with those of the real spacetime, a gauge-invariant description is independent of this identification. Therefore, a gauge- 
invariant approach is conceptionally more attractive although the physical meaning of the gauge-invariant variables 
y \ in terms of which the theory is formulated becomes obvious only for specific observers. 

According to their extension with respect to the Hubble scale, cosmological perturbations are divided into large- 
J^P" scale perturbations, i.e., perturbations with wavelengths larger than the Hubble length, and small-scale perturbations 
^ , for which the opposite is true. 

In inflationary cosmology the presently observed structures are traced back to quantum fluctuations during an 
early de Sitter phase. These originally small-scale perturbations are stretched out tremendously in the inflationary 
period, thereby crossing the Hubble length (which is constant during the de Sitter stage) and becoming large-scale 
perturbations. Lateron, after the inflation has finished, i.e., when the universe is adequately described by the standard 
Fricdmann-Lemaitre-Robertson- Walker (FLRW) model, these perturbations cross the Hubble length again, this time 
inwards, and again become small-scale perturbations. 

Consequently, it is of interest to follow the perturbation dynamics during the time interval between both Hubble 
scale crossings. Bardeen jl]] , Bardeen, Steinhardt and Turner || were the first to recognize the existence of a conserved 
quantity on large scales which simplified the corresponding perturbation analysis. The existence of a conservation 
quantity allows one to establish a link between physical variables in remote cosmological periods. This well-known 
fact has been widely used in the literature (see, e.g., 

Recently the role of conserved quantities was discussed in terms of the gauge-invariant and covariant Ellis-Bruni 
variables || by Dunsby and Bruni [pLl|| . The key quantities of the Ellis-Bruni analysis are the covariantly defined 
spatial gradients of the energy density, the pressure, the expansion and the 3-curvature scalar. From their definition 
all these variables have a general physical meaning for a comoving (with the fluid 4-velocity) observer without any 
reference to perturbation theory. For small deviations from homogeneity and isotropy they may be related to usual 
perturbation variables which are gauge-invariant by construction. The relations among the Ellis-Bruni variables in 
linear order and the corresponding results for the cosmological modes coincide with those of the comoving gauge in 
conventional perturbation theory |T^| . For different gauges, the Ellis-Bruni quantities, if regarded as conventional, 
linear perturbation quantities (see below) have no obvious physical meaning, in general. In other words, as to their 
physical interpretation, comoving observers are preferred. 
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The present paper introduces a set of covariant and gauge- invariant quantities which represent suitable combinations 
of Ellis-Bruni type variables. Like the original Ellis-Bruni quantities, the new variables are general tensorial expressions 
which may be used to characterize deviations from homogeneity and isotropy without explicitly introducing a fictitious 
background universe, i.e., they are no perturbation variables in the conventional sense. While the Ellis-Bruni quantities 
are adapted to comoving hypersurfaces, the new variables are associated to constant expansion hypersurfaces, i.e., 
hypersurfaces of uniform Hubble parameter, hypersurfaces of constant 3-curvature, and constant density hypersurfaces. 

We shall establish the connections of these new quantities both to the Ellis-Bruni variables and among themselves. 
We will reconsider the perturbation theory of a perfect fluid universe and we will show that the corresponding analysis 
formally simplifies if written in terms of the new variables. This simplification is especialy useful with respect to the 
characterization of a conserved quantity in the limit of large perturbation scales. 

For a flat background we find expressions for the latter in terms of the fractional spatial energy density gradient 
(or the spatial gradient of the 3-curvature) on constant expansion hypersurfaces, or, alternatively, in terms of the 
expansion (or 3-curvature) gradients on constant density hypersurfaces. In the case of nonvanishing background 
curvature none of these quantities will be conserved, in general. Only for dust or under the additional condition that 
the curvature terms in the background field equations are small, the above quantities remain approximately constant. 
The most suitable quantity to characterize this case turns out to be the fractional, spatial gradient of the energy 
density on hypersurfaces of constant 3-curvature. 

The paper is organized as follows: In section II we recall the gauge-problem in conventional perturbation theory and 
point out the relation between the gauge-invariant perturbation variables used previously by the author |l4],|l5[ and 
the Ellis-Bruni variables 01. In section III we consider the cosmological dynamics in terms of the Ellis-Bruni variables 
as modified by Jackson fllq ] . In section IV we define covariant and gauge-invariant variables associated with constant 
expansion hypersurfaces, i.e. hypersurfaces of constant Hubble parameter, hypersurfaces of constant curvature and 
constant density hypersurfaces, and we establish relations between them. In section V different expressions for large- 
scale conserved quantities for vanishing background curvature are found and the problems to generalize these results 
to curved backgrounds are discussed. Section VI presents the conclusions of the paper. 



II. THE GAUGE PROBLEM 



The investigation of cosmological perturbations within general relativity was pioneered by Lifshitz | fl7[ . The tradi- 
tional approach (see, e.g. ||) starts with Einstein's field equations 

Gran — Rran ~^QranR — ^mn (f ) 

(m,n.... = 0,1,2,3) and splits both sides of the latter into a zeroth-order background and first-order perturbations 
about this background: 

a =<?(°)-i-r7 t — t(°)+t (i\ 

^ran — ran ' ^ran ; ^ran — - 1 - ran ' ^-ran • V / 

Provided, the zeroth-order problem G™i = nTmn is solved, one has to study the set of perturbation equations 

Gran — ^T mn . (3) 

Restricting ourselves to a perfect fluid with the energy-momentum tensor 

-^mn P^mMn P^ran : (4) 

where p is the energy density, p is the pressure, u m is the 4- velocity and h mn = g mn + u m u n , the equations (||) 
represent a system of differential equations for the perturbation quantities 

P = P~ P {0) , P = P~ P {0) , U m = U m - , g mn = 9mn ~ 9$ n • (5) 

None of these quantities is gauge- invariant, i.e., invariant under infinitesimal coordinate transformations 

x n '=x n -C(x) • (6) 

Assuming a homogeneous and isotropic, comoving (i.e., u m ^ = (1,0, 0, 0),) zeroth order, the transformation properties 
of these quantities are (see, e.g., S, chapter 10.9) 
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p' = p + p^C\ p'=P + P {0) t°, u' m = u m -^ m , (7) 
where we have used p q = p^ = p^iU 11 ^ etc. in the comoving zeroth order, and 

a' =a + C a o (0) + C a o (0) + C a a {0) (8) 

Bran ~ W« ~ S>, raSon ~ S.nilma ^ S iJmn,a • \°) 

One may either chose now a specific gauge, e.g., the traditionally preferred synchronous gauge, or try to find a 
gauge- invariant desciption ( G|). In the latter case one may look for suitable combinations of the quantities (||) that 
are gauge- invariant, i.e., invariant under the infinitesimal coordinate transformations (JgJ) . From the transformation 
properties (Q) it is obvious that, e.g., the quantity 

+ (9) 

(/!... = 1, 2, 3) is invariant under the transformations (^), i.e. pff}' = p$. This quantity represents the energy density 
perturbations on comoving (superscript 'c') hypersurfaces u M = (cf |14|,[l^]). It corresponds to the spatial derivative 
of Bardeen's quantity e m p|Jl3| . 

One easily recognizes that this construction principle of gauge-invariant quantities is applicable for any scalar ( 
|jT4] , fi"5|| ) . E.g., pressure perturbations on comoving hypersurfaces are characterized by 

pf^R.+P^. (10) 

Gauge-invariant descriptions of scalars such as the fluid expansion O = u a a and the spatial 3-curvature 1Z, to be 
defined below, are 

= 0,m + © (0) ^ , ii^ =n^+n^u„ , (ii) 

respectively. All these quantities are gauge-invariant by construction. In order to establish a link between these 
variables and those introduced by Ellis and Bruni Q , we have to investigate the key quantities of the latter authors, 
the spatially projected gradients h c a p c , h c a p c , h c a Q c , h c a lZ. c , etc., in first order. Using the first-order expressions 

h%=0, h° a = u a , (12) 

one finds 

(h c oP J = , {hlp, c ) = + p {0) u, = ffl , (13) 

and corresponding relations for the spatial gradients of p, O, and 1Z in first order. While the quantities h c a p iC etc. 
vanish in zeroth order since the background was assumed to be homogeneous, they coincide in first order with the 
gauge-invariant quantities pf^ etc. 

The second relation ( Jl3| ) demonstrates explicitly how the Ellis-Bruni variables are related to corresponding gauge- 
invariant perturbation quantities with physical interpretations on comoving hypersurfaces. Our considerations may be 
regarded as a kind of motivation for the usefulness of the Ellis-Bruni quantities from the point of view of conventional 
perturbation theory. They also establish the connection beween the latter quantities and the quantities of the type (j^) 
and (|ic|), introduced in |]l^ , |l5| . Using the Ellis-Bruni quantities which, by their definition, vanish in a homogeneous 
universe, one may avoid an explicit decomposition into a fictitious background and perturbations about this (not 
unique) zeroth order. Since the existence of a gauge-problem is just a consequence of the nonuniqueness of the 
background, one circumvents the gauge-problem from the beginning by using the Ellis-Bruni variables which are 
covariant by definition and no perturbation quantities in the usual sense. 

The spatially projected gradients by Ellis and Bruni have a physical meaning on comoving hypersurfaces. As to the 
physical understanding of the gauge- invariant quantities h c a p. c etc., if regarded as conventional first-order perturbation 
variables in the sense described above, the comoving gauge is naturally preferred. In different gauges, these quantities 
have no obvious physical meaning, in general. From these statements the impression may arise that there is something 
special about the comoving gauge, even within a gauge-invariant formalism. This is not the case, however, as we are 
going to clarify in this paper. It is well-known that there exist obviously reasonable gauge-invariant quantities in the 
literature having a simple physical meaning in gauges different from the comoving gauge. For reasons of convenience 
some authors e.g., Bardeen, Steinhardt and Turner ||, occasionally prefer the 'uniform Hubble constant gauge' 
or others, e.g., Hwang [p2| , the constant curvature gauge. One may therefore ask whether there exist covariant 
and gauge-invariant quantities that, from the point of view of their physical interpretation, prefer hypersurfaces of 
constant Hubble parameter or of constant curvature in a similar sense in which the Ellis-Bruni quantities prefer 
comoving hypersurfaces. As we will see below, there are indeed problems for which a choice of covariant variables 
different from the Ellis-Bruni choice may simplify the dynamical description. In order to introduce these variables we 
first recall the basic elements of the cosmological dynamics in terms of the Ellis-Bruni quantities, following here the 
elegant presentation by Jackson J16[]. 
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III. COSMOLOGICAL DYNAMICS IN TERMS OF THE ELLIS-BRUNI-JACKSON VARIABLES 



Instead of applying the field equations (|l]), we shall investigate the cosmological dynamics within the 'fluid-flow' 
approach used in fl^ , jil^ ,fi| pb]j2"l|l . The equations of motion T % \ = 0, imply 

p=-6(p + p) , (14) 

with 

e = , (15) 

and 

(p+p)u m = -p, k h mk , (16) 
where u m = u^u". Additionally, we shall use the Raychaudhuri equation for O, 

9 + Je 2 + 2 (a 2 - lo 2 ) - «» - A + I (p + 3p) = . (17) 
A is the cosmological constant. The magnitudes of shear and vorticity are defined by 

J 2 — ^ ~ ^ ab , , 2 — , , , , ab f1Q\ 

O = -jVabO ) ^ = -jLO ab LO , (1») 



with 



G ab = h c a hfu^ c . d ) - ^Qh ab , uj ab = h a hfu[ c;d ] . (19) 
The 3-curvature scalar of the projected metric, 

K = 2 (-^© 2 -lo 2 + np + A^j , (20) 

reduces to the 3-curvature of the surfaces orhogonal to u a in the case lo = 0. The homogeneous and isotropic FLRW 
universes are characterized by a = lo = ii a = 0. Taking the spatial gradient of (|l4|) yields 

Kp,c = -(p + p) K&, e - QK (p + P ) tC . (21) 

Introducing a length scale S by 

HI (22 » 

and rewriting the l.h.s. of ( pi] ) after multiplying by S as 

Sh c m p, c = h a m (Ship J - eSh c mP . c + « + a c m ) Sh n cP . n , (23) 
the equations ( ^H ) and ( p3| ) may be combined into 

h a n D a + + K + <) D e + t n =0, (24) 

p + p 

where we have introduced the fractional quantity [fl6[ 

_ Sh c a p, c 

D a = ■ , (25) 

p + p 

and 

t a = Sh c a O, c . (26) 
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Equation ( p4| ) is Jackson's eq.(29) (T(|. Similarly, taking the spatially projected gradient of the Raychaudhuri equation 
(|l7|), one obtains 



h a J a = -QP n - K + a c n ) t c - -Qt n - Sh c n (2a 2 - 2uj 2 ) t 
+Sh c n (u° a ) c -~(p+p)[D n + 3P n ] , 



with the fractional quantity 



Sh c a P,c 
P + P 



(27) 



(28) 



characterizing the pressure perturbations. Equation (27) is Jackson's eq.(30) ||lq| . The set of equations ( |24| ) and (|27 
is still completely general. Even for an equation of state p = p (p) which allows one to express P n in terms of D n , the 
equations ( 24) and ( p7|) are, of course, not a closed system for D a and t a since these quantities are coupled to u> and 
a and their spatial gradients. 

From now on we shall assume the spatial gradients as well as a and ui to be small, i.e., we assume the universe 
to be almost homogeneous and isotropic. Consequently, up to first order in the inhomogeneities, the factors in front 
of the quantities D a , P a and t a in (|^) and ( p7j ) refer to the homogeneous and isotropic case with u = a = 0. The 
linearized set of equations becomes 



and 



ll n L a 



KDa 



1 2 

-~np n - -et v 



-D n + t n = , 



P 



where we have used the zcroth-order relations ( A = 0) 

1 



and 



for homogeneous and isotropic universes. The 3-curvature in the latter case is known to be 

n = — , 

with the scale factor a of the Robertson- Walker metric. The last term in ( ^) is generally given by ( fiTi| : 



SK («- ) 



SK 


uab ( una P,m \ 


+ sk 




[ \ a P + PJ; b _ 


,c 



■ % ab P,a Pjb 
p+pp+p 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



In linear order in the inhomogeneities this reduces to 



(35) 



With an equation of state p — p (p) one may write 

p=-cle(p + p) , P n = clD n , (36) 

where c 2 . = (dp/dp) ad is the square of the sound velocity. Because of ( [l3| ) etc., in linear order, the system (|29|), ( |30| ) 
may be written as 



c 2 OD„ 



(37) 
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and 



t„ = -36^ - 



~(p + p) + 3^c? s +c 2 s ^j 



D„ 



(38) 



Eliminating from the system (pTl), (p§), one obtains a closed equation for D^, 



(Cs) ' © + \ (P ~ 3p) C 2 + ~ (p + p) + C 2 S ^ 



r> M = 



(39) 



which corresponds to Jackson's eqation (57). From the latter equation one easily derives the well-known growing 
and decaying modes in the long-wavelength limit for k = 0, equivalent to neglecting the spatial gradient term in the 



bracket in front of D^. 



IV. NEW PERTURBATION VARIABLES 



Let us consider the ratio of the spatial variation of a scalar quantity, say, the energy density, i.e., h c a p^ c to its 
variation in time, p = u c p tC . With (H3) we find, in linear order, 



u c p tC 



P±_ 

P 



(40) 



Obviously, this ratio is a reasonable quantity to characterize small deviations from homogeneity. Similar ratios may 
be formed for the other scalar quantities of interest: 



and 



For k ^ we may also introduce 



U C P,( 



u c e , 



p 



e 



(41) 



(42) 



+ u M {k + 0) 



u c K. r . J K 

For k = 1Z = the 3-curvature perturbations are gauge- invariant, i.e., 



(Kn, 



n, 



(k = 0) 



(43) 



(44) 



By construction, all the quantities (|4(j) - ( |43j ) represent the first-order ratio of the spatial variation of the quantities 
p, p, and 1Z, respectively, to their change in time, on comoving hypcrsurfaces = 0. The 4-velocity perturbation 
enters each of the relations (Eof) 



431) in exactly the same manner, namely simply additively. Forming differences 



between any two of the quantities (^C ) 

([43]) in order to obtain new, covariant and gauge- invariant variables. E.g., we may define 



the velocity perturbations just cancel. This suggests combining (40) 



h c n {ce) 

""mP,c 
u °P,c 



(Jfe = 0,±1) 



(45) 



u P., 



u u p. 



h c m n„ 

U c TZ r 



(k = ±l) 



(46) 



G 



In linear order, with (ffO|), (^) and (^3|), the spatial components are 



' uc ( ce ) \ * A 



^ (fc = 0,±l) , (47) 



and 



u c P,c P @ 



U C P,c P 11 



(k = ±1) . (48) 



Generally, all variables of the type ( f45| ) and ( [46| ) are exact covariant quantities without any reference to perturbation 
theory. As perturbation variables they are gauge-invariant by construction. From ( ff7| ) it seems obvious to conclude 

that ft^p,c / ( uC P,c) in first order represents the ratio of spatial to time variation of the energy density on hypersurfaces 
0. M = 0, i.e., on hypersurfaces of constant expansion or constant Hubble parameter, in the same sense as the Ellis- 
Bruni quantity h c m p, c l (u c p. c ) has a corresponding,, well-known meaning in the comoving gauge = (cf.(|l^ )). 
While we shall indeed use this interpretation of ( |47|) in the following, it requires additional clarification. The point is 
that the quantity p, introduced by (I|), is the energy density for a comoving (with 4- velocity u a ) observer. A different 
observer, moving, e.g., with a 4-velocity n a , normal to hypersurfaces O jjU = 0, would interpret a different quantity, 
namely p = T a \ } n a n b as energy density. (Different from the observer moving with u a he would also measure an energy 
flux.) Generally, p and p are related by (see 

p = p cosh 2 (3 + p sinh 2 f3 , 

where (3(t) is the hyperbolic angle of tilt given by cosh/3 = —u a n a . Corresponding relations hold for the other 
quantities like pressure and expansion. It follows that the perturbation quantity p for a comoving (with u a ) observer 
will generally not coincide with the perturbation p, i.e., p will not coincide with the energy density perturbation on 
hypersurfaces of constant expansion. Now, in the present case we are considering first-order perturbations, while u a 
and n a coincide in zeroth order. For small angles of tilt, however, i.e., for /3<1, the differences between p and p are 
of second order in [3. Consequently, within linear perturbation theory we may identify the quantities p and p. Similar 
statements hold for the other perturbation variables. Therefore, the above interpretation of regarding the quantity 
([l7]) as the ratio of spatial to time change of the energy density on hypersurfaces of constant expansion, M = 0, is 
justified up to first order. Since the latter gauge is defined by the vanishing of the gradient of 0, the quantity 9 itself 
is only determined up to a constant. 

Analogously, for k = ±1, the quantity h^p^? 1 / ( uC p,c) in first order represents the ratio of spatial to time change 
of the energy density on hypersurfaces of constant curvature (superscript 'cc'), defined by 7l tfl = 0. The same 
construction principle may be applied to introduce a variable describing the ratio of spatial to time change of the 
3-curvature, 



u c K. c u c K c u c e . 



(fc = ±I) . (49) 



For k = one has h^TZ^ = hpi., 



The ratio of spatial to time change of the expansion with respect to the corresponding ratio for the 3-curvature is 



u c e c u c e. c u c n , 



(fc = ±l) . (50) 



Again, the 0-component vanishes and the first-order spatial components may be written analogously to ( |47j ) and (^) , 
namely 



and 



(fc = ±l) , (51) 



-di 'A (fc = ±l) . (52) 
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The quantity (51) represents the first-order ratio of spatial to time change of the 3-curvature on constant expansion 
hypersurfaces, while ( |52|) is the corresponding ratio for the expansion on constant curvature hypersurfaces. 

Provided, an equation of state p = p (p) is given, we are left with three quantities being ratios of spatial to time 
changes, namely those for p, and 1Z. Forming differences between them yields variables that in linear order may be 



interpreted as perturbation variables, e.g., on constant expansion hypersurfaces like ( |47| ) and (|5l|), or, like (48) and 
([52|), on hypersurfaces of constant curvature. In a similar way it is possible to introduce variables which in first order 
have physical interpretations on constant density hypersurfaces. The corresponding definitions are 



T->( cd ) h c 1? h c n 



u c 1Z iC u c 7Z, c u c p. c 

and 

h c m Q [c c d) hLQ, c h1p A 



u c Q c u c O jC u c p 
where the superscript 'cd' stands for 'constant density'. In linear order one has 



(k = ±1) , (53) 



(fc = 0,±l), (54) 



= ^ - (k = ±1) , (55) 



and 



u c TZr / n 



K^\ = %E_P^ (fc = o !± l). (56) 



The quantity ( |55| ) is ( f48| ) with the opposite sign. A corresponding relation holds between ( |56| ) and (47). 

The relations of all of these quantities to the Ellis-Bruni variables are obvious. It will turn out that the perturbation 
dynamics looks simpler in any of the new variables than it does in terms of the Ellis-Bruni variables. The reason is 
that the new variables are adapted to hypersurfaces on which one of the first-order spatial gradients of 0, 1Z or p 
vanishes. None of these gradients vanishes, however, on comoving hypersurfaces. 

A. Constant expansion hypersurfaces 

Differentiating the GauB-Codazzi equation ( pp| ) and projecting orthogonal to u a yields, in linear order in the 
inhomogeneities, 

r a = - ^dt a + 2n{p+p)D a , (57) 
where we have introduced the abbreviation 

r a = Sh c a K c . (58) 
According to (^) and ( f47| ) we introduce the fractional quantity 

D^^D a --B-^ (k = 0,±l) , (59) 
P + pQ 

which in linear order represents the fractional, spatial gradient of the energy density on hypersurfaces of constant 
expansion. Analogously, we define the fractional spatial gradient of the curvature on constant expansion hypersurfaces, 

ri C£ ) ee r a 7t| (k = ±1) , (60) 



with (cf @) 



7t=-^9ft. (61) 
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Using the definitions ([59]) and (^) to replace D a and t a , respectively, in (j57|), we find a relation between the gradients 
of the 3-curvature and the fractional energy density on uniform Hubble parameter hypersurfaces: 

rf> =2k( P + p)dM , (fc = 0,±l). (62) 

This is relation (^7|) written in the variables that are adapted to constant expansion hypersurfaces. The latter relation 
is also valid for k — 0, where r a ce ^ ^ r a . 



B. Constant curvature hypersurfaces 

According to (|50| ) and (|52|) the fractional energy density gradient on hypersurfaces of constant curvature is charac- 
terized by 



p r a 

+ pii 



(fc = ±i) 



Analogously, for the gradient of the expansion on these hypersurfaces one has 



e-^ (k = ±i) 

K 



(63) 



(64) 



Inserting D a and t a from ( |63| ) and (|64|) into (|57j ), and using the zeroth-order relations (|3^) and (^Tj), we get 



-e4 cc ) = 2 K ( i0 +p)i?( cc ) (fc = ±i) 



(65) 



which is again (p7|), now in terms of variables adapted to constant curvature hypersurfaces. 

Equation ( |57| ) is a linear relation between the three Ellis-Bruni type perturbation quantities D a , t a and r a defined 
with respect to comoving hypersurfaces. The introduction of the variables ( |59| ) and ( ^KJ) or, (|6^) and (|64|), respectively, 
reduces (l57j ) to relations (either (^) or (|65|)) between only two variables, adapted either to hypersurfaces of constant 
expansion or to those of constant curvature. Effectively, the number of variables has been reduced in both cases. One 
may choose, e.g., one of the fractional energy density perturbations, D a c or Da as independent variable. 

At any stage it is possible to change between the sets of variables associated with either the comoving, or the 
uniform Hubble constant, or the constant curvature hypersurfaces. Subtracting (|63|) and yields 



£)( cc ) _ £)(ce) 



(ce) 

p r a 

p + p ii 



(k = ±i) 



(66) 



Using here the relation (62) between ra and D a ce ' as well as the zeroth-order relations (^||) and (|6ll), we find the 
following connection between the fractional energy density gradients on constant curvature and constant expansion 
hypersurfaces: 



-UD^ = QD a c 



(k = ±1) 



(67) 



C. Constant density hypersurfaces 



We define, according to (^4J) and ( |56| ) with (|26|), the perturbations of the expansion on constant density hypersur- 
faces, 



^t a -(p+p)<d^ (fc = 0,±l), 

P 



as well as, according to (|53|) and (p5|) with (58), curvature perturbations on these hypersurfaces, 



r ^ =r a -{p+p)il^j (k=±l) 



(68) 



(69) 
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Solving eq. <p8\) for t a , eq. (|69[) for r a and inserting into ( |57| ) yields 

r M) = _4@ t ( crf) ( fc = ,±l) (70) 

for the perturbed Gaufi-Codazzi equation in terms of variables that are adapted to constant density hypersurfaces. 
The relation (Foh is also valid for k = 0, where ri cd ' fc =° r a . Introducing into (16^) the quantity by (|59|), one 



obtains a relation between the quantities and £>i° e ^ , 

Analogously, from ( |69| ) and (]6^) one gets 

ri cd) = -\KD^ (k = ±1) . (72) 
It follows, that ta"^ and are connected through 

Qt^ = l -KD^ (k = ±1) . (73) 



V. CONSERVED QUANTITIES 



We are now going to demonstrate that the system (|37j), (|38|) becomes especially simple if rewritten in terms of the 
newly introduced variables. Starting with the definition(|63|) of the fractional energy density gradient on hypersurfaces 
of constant curvature and using @ and dlj), we have 



M cc) = D„ 



3 Ta 

2 K 



(fc = ±l) 



Introducing here (57) for r a and applying the zeroth-order relation (p2[) yields 



£)(<=<=) 



2 r OD a + et 



(fc = ±i) 



(74) 



(75) 



Differentiating the latter relation and using the equations of the system ( |37| ) and (38) on its r.h.s. provides us with 

72 



D 



(cc) _ 



9 -d, - op, - |e^p M (k = ±i) 



p + p 

According to (|3^) the first and second terms on the r.h.s. of ([76] ) cancel and we arrive at 

2 ^V 2 



£)( cc ) 



(76) 



(77) 



This simple relation comprises the entire linear perturbation dynamics for nonflat background universes. In the case 
k — 0, where the 'cc'-quantitics are not defined, one has to modify the above reasoning to arrive at a corresponding 
relation. From the definition (|59|) one finds with (|l3), 



e^i ce) = &D a + et a (k = o, ±1) 



(78) 



The r.h.s. of this equation coincides with the expression in the bracket on the r.h.s. of J7q). After similar steps like 
those between eqs.([75|) and (frj]) one gets 



= -cfQg-z-Dp 
a 2 p 



(k = 0± 1) 



(79) 



Alternatively, taking into account (^32|), (fn]) and (fTOj), the quantity a 2 QD^ may be replaced to yield 
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2 e 



2 K (p + p) 





3 


a a r^ 




~ _ 4 





, 9 V 2 

2 o c 2 — d m 



(80) 



for k = ± 1. For fc ^ equation ([79]) coincides with (|77|), For fc = the quantities a 2 ri ce ' ) and a 2 ri cd ' ) reduce to the 
variable C a used by Dunsby and Bruni jTjJ , 

Equation (^) (or an equivalent formulation according to ( |80| ) or (f77|)) represents the most condensed form of the 
linear perturbation dynamics. It is the main result of this paper. A comparison with ( |39| ) shows that all terms in 
that equation, with the exception of the spatial gradient term, have been included into a first time derivative. The 
price to pay for this formal simplicity is that in general each of the equations (fn]), ([79]) and ( jsc| ) couples variables 
that have physical interpretations in different gauges. 



It is obvious that the introduction of the new variables is especially useful in cases where the r.h.s. of (79) and (JSC 
(or ([Ft] )) may be neglected. 

For = D( n )V nQ( n ) (and corresponding relations for D^ e ^ etc.) where the Qf n ) satisfy the Helmholtz equation 
V 2 Q(„) = — n 2 Q(„) one has ( [pi| , p5P ,pl|l ) n 2 = v 2 for k — and n 2 = v 2 + 1 for k — — 1, where v is continous 
and related to the physical wavelength by A = 2ira/v. For k = +1 the eigenvalue spectrum is discrete, namely 
n 2 = m (m + 2) with m = 1, 2, 3 

On large perturbation scales, i.e., for v <C 1 in a flat background universe the r.h.s. of ( |79| ) and ( |80| ) may be 
neglected and we find 



-(<*) 

a 2 ©^ = a 2 O n P- , 
M 2k(p + P ) 



(81) 



for fc = 0. 



Using equation fl59|) which defines 7?i c e in terms of £> p and and eliminating i M by (|3T|) , yields the following 
relation between and D a : 



-6£) 



■ J ^ 2 



(p + 3p) 



£>( V ) < 1) 



591). On large scales and for k 



(82) 

equation (J82|) together 



Of course, the latter relation and (|79j) are equivalent to 
with is a first order equation to determine Dr v -\ . 

In a nonflat universe (for simplicity we restrict ourselves to the case k — — 1 which is obscrvationally favoured) the 
r.h.s. of (f77f), ( |79| ) and ( p0| ) do not generally vanish even for v <C 1. The only exception is c 2 <C 1, i.e., a dust universe 



(cf |lljl). Using the above eigenvalue structure of the Laplacian in (39) and comparing the last term in the bracket in 
front of in ( ^9| ) with the previous ones by taking into account the zeroth-order equation (|3l|), one finds that the 
spatial gradient term may be neglected for y<l, provided the additional condition 



H 2 a 2 > 1 



(83) 



with H = 0/3 is satisfied. The latter relation, however, coincides with the condition under which the curvature 



term in the background equations ( |3l| ) and (|32|) may be neglected. Consequently, the quantities D^ c \ a 2 ©!? 
and equivalent expressions under the dot derivative on the l.h.s. of (^0|) are only conserved on large scales if the 
background curvature terms are negligible. Corresponding properties for a differently defined quantity were found 
by Dunsby and Bruni [jll[(see eq(16) in pT| . For special cases these authors, however, also constructed conservation 
quantities without the restriction ( p3| ) and on arbitrary scales.) 

While the quantities D^ e \ r^ e \ t^ d ^ and r^ d ^ which appear on the l.h.s. of ( f79| ) and ( [so] ) are well defined for any 

(cc) 

value of the background curvature, the quantity , the fractional energy density perturbation on hypersurfaces 
of constant curvature, defined only for k ^ 0, appears to be the most suitable quantity to characterize perturbations 
on a nonflat background, however small the curvature terms may be. The general perturbation dynamics reduces to 
((77]) in this case and is conserved on large scales if a 2 H 2 3> 1. 



ce) 



VI. CONCLUSIONS 



We have simplified the linear, perfect fluid cosmological perturbation theory by the introduction of new, covariant 
and gauge-invariant perturbation variables. These variables are physically interpreted on hypersurfaces of constant 
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expansion, constant 3-curvature, or constant energy density. Their relations to the Ellis-Bruni variables which have 
physical interpretations on comoving hypersurfaces are established. The conserved quantity in the large-scale limit for 
k = may be expressed in terms of the gradient of the fractional energy density on constant expansion hypersurfaces 
or, alternatively, in terms of the expansion gradient on constant density hypersurfaces. The most suitable quantity for 
the description of the perturbation dynamics in a universe with nonvanishing background curvature is the covariantly 
defined fractional, spatial gradient of the energy density on hypersurfaces of constant 3-curvature. It is this quantity 
which is approximately conserved on large scales provided the background curvature terms are small. 
One may hope that use of the variables introduced in the present paper also simplifies perturbation calculations under 
more general circumstances. 

Acknowledgment 

This paper was supported by the Deutsche Forschungsgemeinschaft. I thank both referees for helpful suggestions and 
comments. Discussions with Diego Pavon and Jai-chan Hwang are gratefully acknowledged. 



Bardeen J M 1980 Phys.Rev. D22 1882 

Kodama H and Sasaki M 1984 Prog.Theor.Phys.Suppl. 78 1 

Ellis G F R and Bruni M 1989 Phys.Rev. D40 1804 

Stewart J M 1990 Class. Quantum Grav. 7 1169 

Weinberg S 1972 Gravitation and Cosmology (New York: Wiley) 

Bardeen J M, Steinhardt P and Turner M 1983 Phys.Rev. D28 679 

Lyth D 1985 Phys.Rev. D31 1792 

Hwang J 1991 Astrophys.J. 380 307 

Mukhanov V F, Feldman H A and Brandenberger R H 1992 Phys.Rep. 215 203 

Liddle A R and Lyth H 1 1993 Phys.Rep. 231 1 

Dunsby P K S and Bruni M 1994 Int.J.Mod.Phys. D3 443 

Ellis G F R, Hwang J and Bruni M 1989 Phys.Rev. D40 1819 

Bruni M, Dunsby PKS and Ellis G F R 1992 Astrophys.J. 395 34 

Zimdahl W 1991 Class. Quantum Grav. 8 677 

Zimdahl W 1992 Astrophys.J. 393 471 

Jackson J C 1993 Mon.Not.R.Astron.Soc. 264 729 

Lifshitz E M 1946 J. Phys. (USSR) 10 116 

Hawking S W 1966 Astrophys.J. 145 544 

Olson D W 1976 Phys.Rev. 14 327 

Lyth D H and Mukherjee M 1988 Phys.Rev. 38 485 

Lyth D H and Stewart E D 1990 Astrophys.J. 361 343 

Hwang J 1993 Phys.Rev. D48 3544 

King A R and Ellis G F R 1973 Comm. Math. Phys. 31 209 
Lifshitz E M and Khalatnikov I M 1963 Adv. Phys. 12 185 
Harrison E R 1967 Rev.Mod.Phys. 39 862 



12 



